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ABSTRACT 

The spatial variation of physical quantities, such as the mass density, across solar atmospheric waveguides 
governs the timescales and spatial scales for wave damping and energy dissipation. The direct measurement 
of the spatial distribution of density, however, is difficult and indirect seismology inversion methods have been 
suggested as an alternative. We applied Bayesian inference, model comparison, and model-averaging tech¬ 
niques to the inference of the cross-field density structuring in solar magnetic waveguides using information on 
periods and damping times for resonantly damped magnetohydrodynamic (MHD) transverse kink oscillations. 
Three commonly employed alternative profiles were used to model the variation of the mass density across the 
waveguide boundary. Parameter inference enabled us to obtain information on physical quantities such as the 
Alfven travel time, the density contrast, and the transverse inhomogeneity length scale. The inference results 
from alternative density models were compared and their differences quantified. Then, the relative plausibility 
of the considered models was assessed by performing model comparison. Our results indicate that the evi¬ 
dence in favor of any of the three models is minimal, unless the oscillations are strongly damped. In such a 
circumstance, the application of model-averaging techniques enables the computation of an evidence-weighted 
inference that takes into account the plausibility of each model in the calculation of a combined inversion for 
the unknown physical parameters. 
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1. INTRODUCTION 

The dissipation of magnetic wave energy is considered a 
relevan t process in the heating of the plasma in the solar atmo¬ 
sphere (iParnell & De MoorteJ2012tlArreguil2015h . High res¬ 
olution observations show tha t wave activity is all pervasive in 
the solar atmosphere (see e.g.| Nakariako^& yer wicht3l2005l : 
De Moortel & Nakariakov^l2fljess et alJl2015h . but quanti- 
fying the importance of magnetohydrodynamic (MHD) waves 
in plasma heating processes remains a difficult task. The sig¬ 
nificance of currently considered physical processes leading 
to wave energy transport and dissipation s uch as resonant ab¬ 
sorpti on (lGoossenslll991h . phase mixi ng (iHevvaerts & Priest 
1^^), or Alfven wave turbulence (Ivan Ballegooiien et al. 
201 ll) relies on the variation of physical properties across 
the magnetic field. In particular, the cross-field variation of 
density and magnetic field strength determines the timescales 
and spatial scales for the damping of tran sverse waves 
(iHollweg & Yangjll988l: iGoossens et al.ll201 Ih. how fast en¬ 

ergy is transferred to srnall len gth scales (iHeyyaerts & Priesd 
119831: ISoler & TerradasI 1201 5b. the timin g for the onset of 
dissinatiye effects (iLee & Robertsl [Hl^, and the fraction 
of the waye energy ayailable to be conyerted into heat 
(IGoossens et alJl2013l : IVan Doorsselaere et alJl2014!) . 

Seismology inyersion techniques offer an indirect method 
to obtain information on the physical properties of mag¬ 
netic and plasma structures and their spatial yariati on. The 
method of MHD seismology, first suggested by lUchidal 
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(119701) . iRosenbei^ (Il970l) . and iRoberts et al.l (Il984h . com¬ 
bines information on obseryed and theoretical properties of 
MHD wayes to obtain information on difficult-to-measure 
physical parameters. Early implementations of the method 
haye proyen to be useful in the determination of e.g., the mag- 
netic field strength in transy e rsely oscillating coronal loop s 
(iNakariakoy & OfmanI |200H IVan Doo rsselaere et al.l 1200^ . 
the coronal density scale height dAndries et al j 2^00.51) . the 
magnetic flux tube expansion (iVerth et aU 2008 |). or th e 
Alfyen yelocity i n pro m inence threads ( On^^l] 120091) . 


see iDe Moortell __ 

iDe Moortel & Nakariako' 


(l200l. In 

cariakoy| (l2 


(Nakariakoy & Verwichtd (120051) . 
2 OI 2 I) , an d lAireguinEoiTh for re- 


yiews. 

A particular application of MHD seismology aims at ob¬ 
taining information on the spatial yariation of the mass den¬ 
sity across solar coronal wayeguid es. An early example can 
be found in IGoossens et al.l (l2002h and inyolyes the compari¬ 
son between the theoretically predicted damping ratio of res¬ 
onantly damped transyerse kink wayes in radially inhomoge¬ 
neous wayeguides and the obseryed damping ratio of coro¬ 
nal loop oscillations. By assuming a density contrast be¬ 
tween the interior of the loop and the external corona, esti¬ 
mates for the transyerse inhomogeneity length scale can be 
obtained. Recent expansions of this idea haye considered 
the Bayesian inyersion of both density contrast and trans¬ 
yerse inhomogeneity l ength scale by mak ing use of differ¬ 
ent damping regimes (lArregui et al.ll2013h or the computa¬ 
tion of their marginal posteriors from obseryed damping ra¬ 
tios (lArregui & Asensio Ramos|l2014l) . In these studies, the 
density structuring was ar bitrarily pre s cribed . 

In t wo recent papers, ISoler et ^ (120131) and ISoler et ^ 
(120141) haye considered the impact of the assumed cross-field 
density profile on theoretically predicted periods and damping 
times for resonantly damped transyerse oscillations and the 
ensuing in fluence on the inferr ed physical parameters through 
inyersion. ISoler et aH ( 120131) describe the basic theory be- 
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hind the application of the Frobenius method implementa¬ 
tion to compute resonantly damped oscillat ions in nonuniform 
waveguides. The seismology analysis by ISoler et aH (1201 4l) 
for parameters such as the Alfven speed, the density contrast, 
and the transverse inhomogeneity length scale shows that sig¬ 
nificantly different inversion results are obtained depending 
on the adopted density model. 

In our study, a comparative analysis be t ween t he alternative 
density models employed bv ISoler et al.l (1201 4h is presented. 
The aim is to quantify to what extent they result in distinct 
parameter inference results and to evaluate which one among 
the presented models better explains observed data for period 
and damping times. This is done by considering a Bayesian 
approach to the solution of the inference problem and by com¬ 
puting the relative plausibility between models. We addition¬ 
ally present a model-averaging procedure to compute a com¬ 
bined inversion result, weighted with the evidence for each 
model in view of observed data. 

The layout of the paper is as follows. Section |2] describes 
the models considered in our study of transverse oscillations 
in one-dimensional flux tubes with alternative forms for the 
cross-field variation of the density. In Section[3]the solutions 
to the forward problem for transverse oscillations and the clas¬ 
sic inversion results are discussed. Parameter inference results 
in the Bayesian framework are described in SectionlH Then, 
the plausibility between the considered alternative models is 
discussed in Section|5] Evidence-weighted posterior distribu¬ 
tions for the unknown parameters are obtained in Section |6l 
In Section|2l our conclusions are presented. 

2. CROSS-FIELD DENSITY MODELS 



r/R 


Figure 1. Spatial variation of the equilibrium mass density across the mag¬ 
netic flux tube for waveguide models with a linear profile (solid line), a sinu¬ 
soidal profile (dotted line), and a parabolic profile (dashed line). The ratio of 
internal to external density is pi/pe = 10 and the thickness of the nonuniform 
layer is IjR = 1. 

density modeR _ 

Following ISoler et alJ ([2013), we consider three alternative 
models for the density variation across the magnetic field. 
We use , M^, and to name the sinusoidal, linear, and 
parabolic density models used in this work. By adopting the 
same superscript convention for the radial variation of the 
density at the nonuniform layer, we have 



The simplest models that still capture global properties of 
MHD kink waves in solar flux tubes, such as the oscillatory 
period and the damping time by resonant absorption, con¬ 
sider the waveguides as one-dimensional density enhance¬ 
ments in cylindrically symmetric models with the magnetic 
held pointing along the z-direction. Under the zero plasma-/! 
approximation, the density distribution is arbitrarily imposed 
to model the waveguide as a region of enhanced mass den¬ 
sity that varies in the radial direction to connect the internal 
higher-density medium to the relatively lower external den¬ 
sity region. The transition between the internal, pi, and the 
external, p^, densities occurs at a nonuniform layer of thick¬ 
ness IfR, with R the radius of the waveguide. At a given 
location inside the nonuniform density layer the global kink 
mode frequency matches the local Alfven frequency. This 
produces the ti me decay of the wave amplitude for stand¬ 
ing kink waves (iHollweg & Yangll988tlGoossens et al.ll2002t 
Ruderman & RobertsI |20021 l This resonant damping pro¬ 

cess also operates for propagating kink waves, resulting in 
the spatial decay of the wave amplitu de as the wave prop¬ 
agates guided by the plasma structure (iTerradas et al.ll20Tot 
iPascoe et al.ll2012h . 

Initial theoretical investigations on the damping of trans¬ 
verse kink waves by resonant absorption focused on the con¬ 
sideration of thin tube and thin boundary models. This is the 
so-called TTTB approximation, which assumes that the wave¬ 
length is long and that the thickness of the nonuniform transi¬ 
tional layer is short in comparison to the radius of the waveg¬ 
uide {rJa <g 1, IjR 1) ('see, e.g., iHol l weg & Yang|[ 19881 
iGoossens et al.|[T995Ll200^jRuderman & Robertsll2002l) . The 

analysis was then extended to include fully nonuniform lay¬ 
ers in both one- dimensional (IVan Doorsselaere^UaL 12004 
I Soler et ^120131) and two-dimensional dArregui et al 1200^ 


for the sinusoidal model. 


pV)=Pi-^^(T-R+^j, (2) 

for the linear model, and 

P^(r) = Pi - 

for the parabolic model. None of the specifically assumed 
profiles is expected to be an accurate quantitative represen¬ 
tation of the real density variation across solar flux tubes, 
but they provide a means to model alternative cases that, by 
means of the model comparison technique presented in this 
paper, can offer information about the plausibility of different 
profiles. Figure [1] shows the spatial variation of the density 
distribution for the three considered profiles for fixed values 
of the density contrast and the thickness of the nonuniform 
layer. 


3. FORWARD AND INVERSE SOLUTIONS 

The basic theory for the time/spatial damping of stand¬ 
ing/propagating kink waves due to resona nt absorption has 
been d eveloped by a numbe r of studies (IHollweg & Yang! 
1988t IGoossens et alJ 1 19951 1201 ll [Terradas et alJ 1201 Ot 

alc~ 


Verth et al.l 120 lOt ISoler et ^ 1201 lald blL Numerical 


ulations have confirmed the obtained damping properties, 
by analyzing the spatial and temporal characteristics of the 
mode coupling process in corona l loop models and a r¬ 
bitrary inhomogeneous structures (ITerradas et alJ l2006allR 
l2008tlPascoe et alJ2010ll201 lLl2012h . For standing waves, an 
analytical expression for the period of the fundamental kink 
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mod e can be obtained und er the thin tube approximation (see, 
e.g. JGoossens et an2008h 

= (4) 

Here ta.i is the internal Alfven travel time and ^ = Pi/Pe 
the density contrast. This solution to the forward problem 
expresses a relation between the observable quantity (pe¬ 
riod) and two physical parameters to be determined (Alfven 
travel time and density contrast). The period is indepen¬ 
dent of the transverse inhomogeneity length scale, IjR, be¬ 
cause of the adopted thin tube approximation. It becomes de- 
pendent on this parameter when this assumption is relaxed 
(IVan Doorsselaere et al.ll200^lArregui et al.ll2005h . The error 
on the period associat ed with the use of the TTTB approxi¬ 
mation was studied by ISoler et ^ (120141) 

Because of resonant absorption, time damping occurs and 
the amplitude of the radial velocity component decays, while 
the azimuthal velocity component becomes larger in the vicin¬ 
ity of the resonant position. The classic theory for reso¬ 
nant Alfven waves under the thin boundary approximation 
provides us with an analytical expression for the damping 
time, Td, as a function of the relevant physical parameters. 
In units of the oscillatory period, thi s expression is of the 
form (see, e.g.,[Hollweg & Yang|ll988l:lGoossens et d] |200l 

Ruderman & Robertsll2002l) 



The numerical factor F depends on the assumed density pro¬ 
file at the nonuniform layer. The factor is Ijn for the sinu¬ 
soidal profile, 4/;7r^ for the linear profile, and 4 V2/7r^ for the 
parabolic profile. The numerical values for these factors do 
not seem to differ to a great extent, but the analytical inver¬ 
sions performed bv ISoler et al.l (1201 4l) indicate that significant 
differences may arise between the resulting one-dimensional 
inversion curves when fully numerical results to the forward 
problem are used. It must be noted that analogous expres¬ 
sions to the ones given in Eqs. (|4|i and d^i can be obtained 
in the case of spatia l damping of propagating kink waves 
(iTerradas et al.l[2Mol) . upon replacement of the period by the 
wavelength and the damping time by the damping length. The 
reason is that resonant damping does not make any distinc¬ 
tion with respect to the standing or propagating character of 
the wave and, under the thin tube and thin boundary approx¬ 
imations, the influence of the particular density profile at the 
nonuniform transitional layer is solely contained in the fac¬ 
tor F. In bot h cases, the dynamics corres ponds to a surface 
Alfven wave (iGoossens et al.l2()0^l2012ah . For standing kink 
waves, the observational consequence is the attenuation of the 
amplitude in time. For propagating kink waves, resonant ab¬ 
sorption produces the attenuation of wave amplitude in space. 

Once the solution to the forward problem is obtained by 
solving Fqs. (|4]i and (|5ll for a range of values of the phys¬ 
ical parameters, their inference can be perfo rmed using the 
analyt ical inversion technique developed by IGoossens et alJ 
(l2008l) . In the inversion procedure, the three quantities (tai, 
and IjR) are the seismic quantities of interest to be inferred 
from observed period and damping time. As we only have 
two observables and three unknowns, there are in principle 
an infinite number of solutions that ca n equally well explai n 
the observed data, as first pointed out bv lArregui et al.l (120071) . 



Figure 2. Classic seismological inversion result of the event with strong 
damping (P = 185 s and = 200 s) with the one-dimensional solution 
curves in the three-dimensional parameter space corresponding to the sinu¬ 
soidal (black), linear (red), and parabolic (blue) density models. The dashed 
lines are the projections of the inversion curves to the various planes. 

These solutions follow a one-dimensional curve in the three- 
dimensional parameter spac e of unknowns. Th e curve was 
first obta ined numerically bvlArregui et aP (120071) and analyt¬ 
ically bv iGoossens et al.l (120081) for a particula r density model 
in the c ontext of standing kink waves and by IGoossens et al.l 
(I2012bl) in the context of propagating kink waves. The 
method is applicable to both coronal loop and prominence 
fine-structure oscillations (see rev iews by iGoossenj 120081: 
lArregui et al.l 1201^ lArregMl2012l) . Solution curves for the 
three densi t y mod els under consideration were obtained by 
ISoler et ^ ( 120141) . Although the inversion solution is unable 
to fully constrain the three parameters of interest, the seis¬ 
mic variables are found to be constrained to intervals that 
can be calculated analytically, under the thin tube and thin 
boundary approximati ons. Their expression s can be found in 
IGoossens et al.l (120081) and iGoossensI ([200^. 

The forward solutions for period and damping time, out¬ 
side the thin tube and thin boundary approximations, have 
been computed for par t icular cross-field density models by 
IVan Doorsselaere et al.l (l2004l) in one-d imensional equilib- 
rium models and bv lArregui et al.l (120051) in two-dimensional 
equilibrium models. The two observables are the n functions 
of the three unknown parameters, tai, and IjR. ISoler et Tl 
(120141) computed numerical solutions for the three density 
models described above and compared the differences be¬ 
tween he numerical forward solutions and those computed un¬ 
der the TTTB approximations and also between the forward 
solutions obt ained with the thre e alternative density models. 
Furthermore, ISoler et ^ (120141) applied the analytical inver¬ 
sion procedure to forward solutions obtained using the three 
density models. 

An example inversion solution (see Figure |2|l, shows the 
one-dimensional solution curve that links those valid combi¬ 
nations of (tai, f, UR) that reproduce a particular observed 
period and damping time. As mentioned in that study, the im¬ 
pact of the assumed density profile at the nonuniform layer 
is strong both in the shape of the one-dimensional inversion 
solution curves in the three-dimensional parameter space and 
on the obtained valid intervals of the seismic variables. Thus, 
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the inferred intervals of the seismic variables are directly af¬ 
fected by the sp ecific choice of de nsity variation. From the 
obtained results. [Soler et al.) (l2014l) conclude that the particu¬ 
lar choice for the density model at the nonuniform layer can 
have a strong impact on the seismologically inferi'ed parame¬ 
ters. 

In our study, we will employ both analytical approxima¬ 
tions (given by E qs. PI and |[5l) an d the numerical solu¬ 
tions, as given by I Soler et ^ (l2014h . to the forward prob¬ 
lem. The inversion procedure, however, will differ, and, in 
addition, model comparison and averaging techniques are im¬ 
plemented. The classic approach to the inversion problem 
deals with finding the mathematical solution for the one¬ 
dimensional inversion curve in the three-dimensional param¬ 
eter space by imposing the exact matching of the theoreti¬ 
cally predicted period and damping time values with those 
observed. The procedure being mathematically correct, we 
recall that parameter inference has to be pursued under cir¬ 
cumstances in which information from observations is incom¬ 
plete and uncertain. For this reason, the inference has to be 
cari'ied out by following a probabilistic approach. We will 
apply three levels of Bayesian inference to the problem of de¬ 
termining the cross-field density profile from observed values 
of period and damping time. The first level involves the in¬ 
version of the three parameters of interest, the internal Alfven 
travel time, and the two parameters that define the cross-field 
density profile, namely, the density contrast and the transverse 
inhomogeneity length scale, assuming that each of the mod¬ 
els described above is true. The second level deals with the 
comparison between the three assumed models, to ascertain 
which one is more plausible in view of given data. The third 
level uses all the information on the data and the three models 
under consideration to produce parameter inversions in which 
the resulting inference is an average of the particular infer¬ 
ences for each model, weighted with the plausibility for each 
model in view of data. 


vides us with the probability density function of the unknown 
parameters conditional on the observed data and the model 
M, p{0\d,M), as a combination of the likelihood function, 
p{d\0, M), and the prior probability density function for the 
unknown parameters, p(8\M). Their relationship is given by 
the expression 


p{0\d, M) 


p{d\e, M)p{e\M) 
f p(dl0,M)p(0lM)d0’ 


( 6 ) 


which assumes that model M under consideration is true. The 
denominator is the so-called evidence, an integral of the likeli¬ 
hood over the prior distribution that normalizes the likelihood 
and turns it into a probability. This can be ignored when per¬ 
forming parameter inference, because it is independent of the 
vector of parameters, but plays the central role in model com¬ 
parison, as will become apparent in Section|5] Both prior and 
likelihood represent probabilities that are directly assigned, 
whilst the posterior is computed. This posterior contains all 
the information that can be gathered concerning the vector of 
unknowns, conditional on observed data. The magnitude of 
the posterior probability density function is a measure of the 
degree of belief on the values the parameter vector can take 
on. 

We proceed by first assigning the direct probabilities to be 
used in Bayes’ theorem for the likelihood function and the 
prior distribution. We adopt a Gaussian likelihood function 
that relates the observed period and damping time, (F, tj), 
and the predictions of the model, (P^(6), t^(6)). Under the 
TTTB solutions they are given in Eqs. (|4| and Q, where the 
particular value of F corresponding to each density model has 
to be used. Outside the TTTB approximation they are nu¬ 
merically computed. Then, the observed values and the theo¬ 
retical/numerical predictions can be compared by adopting a 
likelihood of the form 


4. PARAMETER INFERENCE 

Eigure |2] indicates that the classic inversion leads to well- 
differentiated results for the inference of the three unknowns 
depending on the assumed cross-field density model. Our 
first aim is to evaluate in a quantitative manner how different 
the results are when a fully Bayesian approach is employed, 
which, in addition, enables us to coiTectly propagate uncer¬ 
tainty from measured wave properties to inferred parameters. 

We have performed the inference for the unknown phys¬ 
ical parameters, tai, and l/R by first generating theoreti¬ 
cal predictions for period and damping time by resonant ab¬ 
sorption. The latter two parameters completely define the 
cross-sectional density profile, once one of the three alter¬ 
native models (linear, parabolic, sinusoidal) is selected. The 
forward predictions by each model are computed for different 
combinations of the equilibrium parameters both under the 
TTTB approximations, using the algebraic expressions © 
an d dSll, and by emp loying the Erobenius method developed 
by ISoler et ^ (1201 3h and applicable to any arbitrary density 
profile at the nonuniform layer. In the following, we will refer 
to them as the TTTB and the numerical forward solutions, re¬ 
spectively. Those predictions are then compared to observed 
data for period and damping time in the inversion. 

Let us gather the three parameters to be inferred in the 
vector of unknowns 0 - (tai,^, and the two observ¬ 
able quantities in = (P , Td). The inference is based on 
the use of Bayes’ theorem (iBaves & PricellT76^ . which pro¬ 


piP,Td\0,M) = {Ina-po-rJ ' 

P _ 

2cr], 

with cTp and ctt^ the uncertainties associated with the mea¬ 
sured period and damping time, respectively. Concerning the 
prior information, we adopt uniform prior distributions for the 
unknown parameters over given ranges, so that all the values 
inside those ranges are equally probable a priori. We therefore 
assign 


xexp< 


( 0 )]' [rd-r^{e)f 


2<t^ 


■a) 


p(0i) = 


1 


^ax _ 


for 0”™ < e,-< 0”’“ i = 1,2,3 


( 8 ) 

and zero otherwise, where 9\, 02 , and 63 correspond to tai, 
and l/R, respectively. Unless otherwise stated, the following 
ranges have been taken: tai e [1,500] s, ^ e [1.1,10], and 
l/R e [0.01,2.0]. Application of Bayes’ rule to expressions 
(|7| and ([^ provides us with the full posterior, p(0\d, M). 

In our analysis, rather than sampling the posterior 
using a Markov Chain Monte Carlo technique as in 
Arregui & Asensio RamosI (1201 ll) . we make use of the ba¬ 
sic definition of marginal posteriors, by performing the re- 
qu ired integrals over the parameter s pace. This was first done 
bv l Arregui & Asensio RamosI (l2014l) to infer the density con- 
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trast and transverse inhomogeneity length scale from mea¬ 
sured damping ratio values. As shown by these authors, al¬ 
though the classic one-dimensional curves in the parameter 
space point to an infinite number of equally valid solutions, 
not all of them are equally probable. To obtain the degree of 
belief in terms of marginal posteriors, the following integrals 
have to be computed; 

piTAi\{P, Td), M) = j p({TAi, l/R]\{P, Td), M) d(llR), 

p((\{P, Td), M)^ JpUtai, l/R]\{P, Td), M) dTAi d(l/R), 

and 

pil/R\{P, Td), Af) = Jp({TAi, l/R]\{P, Td), M) dTAi d(- 

(9) 

The marginal posteriors encode all the information for each 
model parameter available in the priors and the data. Their 
computation correctly propagates uncertainty from data to in¬ 
ferred parameters. The a nalysis is applied to two cases from 
a set of observations by lOfman & AschwandenI (l2002h . one 
representing observations with moderate damping (P - 272 
s, Td = 849 s, Td/P = 3.12) and another representing ob¬ 
servations with strong damping {P - 185 s, Td = 200 s, 
Td/P = 1.08). A timescale of 30 s is taken as the uncertainty 
on measured period and damping time, in line with current 
observational capabilities. 

For the case with moderate damping, application of Bayes’ 
rule, using the likelihood function in Eq. 0, uniform prior 
distributions, and the marginalization rules given in Eqs. 0 
leads to the marginal posterior density functions displayed in 
Eigs.[3]and|4] In Eig. [3 the results are arranged in such a 
way that a comparison between the inversion performed for 
TTTB results and fully numerical results can be directly com¬ 
pared. Each of the three columns displays results for each of 
the three used density models. The Alfven travel time can be 
properly constrained (see Eigs. [3^-c). There are no signifi¬ 
cant differences depending on whether we employ the analyt¬ 
ical approximations or the full numerical forward solutions. 
The inversion suggests that low density contrast values are 
preferred over large contrast ones (see Eigs. |3}l-f)- However, 
the posteriors for density contrast display a rather long tail. 
This leads to this parameter being constrained with a large un¬ 
certainty. More constrained marginal posteriors are obtained 
for the transverse inhomogeneity length scale (see Eigs. [3^- 
i). They show that short values of IjR, below IjR - 0.5, are 
more plausible than models with a fully nonuniform layer. 
The largest discrepancies between the TTTB approximation 
results and the fully numerical results are obtained for the in¬ 
ference of the density contrast made by employing the linear 
and parabolic profiles. Eor both the Alfven travel time and the 
transverse inhomogeneity length scale, analytical and numer¬ 
ical forward solutions led to similar inversion results. 

In Eig.m the results are arranged in such a way that a com¬ 
parison between the inversions for the three alternative den¬ 
sity models can be directly made, with the top row showing 
TTTB results and the bottom row numerical results. Eor the 
inference of the Alfven travel time, using the TTTB approxi¬ 
mations for the forward problem, the inference is almost iden¬ 
tical regardless of the employed density model (see Eig. |4^). 
The density model only affects the forward solutions through 
the factor F in Eq. 0; hence the impact is minimal on the in¬ 
ference of the Alfven travel time and only due to the coupling 


with Eq. 0 through the density contrast parameter. As the 
damping is moderate, slightly different marginal posteriors 
are obtained when fully numerical solutions for the forward 
problem are employed (see Eig.lHl). Eor the density contrast, 
outside the TTTB approximation, changing the density model 
slightly affects the shape and probability density functions for 
different density models (see Eig.0). The largest differences 
in the computed marginal posteriors are obtained for the trans¬ 
verse inhomogeneity length scale (see Eigs. |4}; and f). Eor 
this parameter, differences in the marginal posteriors depend¬ 
ing on the density model already appear when the analytical 
TTTB approximations are employed. 

Table [3 summarizes the marginal posterior density func¬ 
tions obtained from the inference for the three density mod¬ 
els and the analytical and numerical forward solutions. Eor 
each parameter and case, the median of the marginal proba¬ 
bility density function and errors at the 68% credible inter¬ 
val are given. Despite the differences in the results discussed 
above, the results indicate that when the posteriors are sum¬ 
marized, almost the same results are obtained, regardless of 
the method employed to solve the forward problem and the 
adopted density model. We can therefore conclude that, for 
moderate damping, the use of either TTTB or numerical for¬ 
ward solutions leads to very similar inversion results and that 
the impact of the assumed density model is small. 

The analysis is now repeated for the case with strong damp¬ 
ing. The corresponding marginal posteriors are displayed in 
Eigs. |5] and |6] In Eig. |5] we compare the results between the 
inversion performed using TTTB and numerical forward so¬ 
lutions. The inference for the Alfven travel time is similar in 
both cases, when a sinusoidal density model is employed (see 
Fig-0), but differs when either the linear or the parabolic den¬ 
sity models are employed (see Eigs. 0 and c). In these two 
cases, the inference with the TTTB forward solutions leads to 
probability density functions that are shifted toward smaller 
values of the parameter. The inference for the density con¬ 
trast leads to similar results between the TTTB and numerical 
cases, when the sinusoidal density model is used (see Eig.0), 
while it shows more marked differences when either of the 
other two profiles is adopted (see Eigs.0 and f). The largest 
differences between the inference applied to TTTB and nu¬ 
merical forward solutions is obtained for the transverse inho¬ 
mogeneity length scale (see Eigs.0-i), but in this case they 
are less important in the case of the parabolic density model. 

In Eig. |6] the results are arranged in such a way that 
a comparison between the inversions for the three alterna¬ 
tive density models can be directly made, with the top row 
showing TTTB results and the bottom row numerical results. 
Eor the inference of the Alfven travel time and using the 
TTTB forward solutions the same probability density function 
is obtained, regardless of the employed density model (see 
Eig. 0). Slightly shifted marginal posteriors are obtained, 
when numerical forward solutions are employed (see Eig.|6}l). 
We note here that the period is also a function of the thick¬ 
ness of the layer when no analytical approximations are used. 
The shifted marginal posteriors for the Alfven travel time in¬ 
dicate that the impact of the density model is measurable in 
the inference of this parameter. Eor the density contrast, the 
effect of the adopted density model is also unimportant under 
the TTTB approximations (see Eig. 0), but more significant 
differences are visible between the marginal posteriors ob¬ 
tained from the numerical forward solutions (see Eig.0). If 
we compare the inversions for the density contrast performed 
from numerical solutions for the moderate and strong damp- 
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Figures. Marginal posteriors for the three parameters of interest (tai, l/R) computed using expressions for the moderate damping case {P = 272 s; 

= 849 s; (Tp = (Tr^ = 30 s). Left column: sinusoidal density profile; middle: linear profile; right: parabolic profile. In each panel, solid lines represent the 
inference performed using the TTTB forward solutions, while dashed lines are for the inference performed using the numerical forward solutions. Summary data 
for the posteriors with the median and errors at the 68% credible region are shown in Tabled 
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Figure 4. Marginal posteriors shown in Fig.l^now displayed for direct comparison between results for the three alternative density models: sinusoidal (solid), 
linear (dotted), parabolic (dashed). Top row: TTTB results; bottom row: numerical results. 
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Figure 5. Marginal posteriors for the three parameters of interest (tai, IjR) computed using expressions m for the strong damping case (P = 185 s; ra = 200 
s; (Tp = (Tr^ = 30 s). Left column: sinusoidal density profile; middle: linear profile; right: parabolic profile. On each panel, solid lines represent the inference 
performed using the TTTB forward solutions, while dashed lines are for the inference performed using the numerical forward solutions. Summary data for the 
posteriors with the median and errors at the 68% credible region are shown in Table[2] 
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Figure 6. Marginal posteriors shown in Fig.[5]now displayed for direct comparison between results for the three alternative density models: sinusoidal (solid), 
linear (dotted), parabolic (dashed). Top row: TTTB results; bottom row: numerical results. 
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Table 1 

Summary of inference results for the moderate damping case for the 
three alternative models by employing the TTTB and numerical 
forward solutions. P = 272 s; Td = 849 s; crp = = 30 s. 


s, 

TTTB approximation 


Numerical 



M^ 

M^. 

mP 

M® 

M^ 

Mf 


lJl +22 

^ ‘ -24 

169+24 

171+21 
^ ' -24 

177+24 
^ ^ -20 

169+2« 



4 4 + 3.7 
- 2.6 

4 1+3-9 

H-.i-2_2 

4 3+3.8 
^•^-2.6 

7 q + 8.2 
'•^- 6.1 

4 . 6 +^°/ 

6-7^8 

l/R 

^•■^- 0.1 

- 0.1 

^•■^- 0.1 

^•^- 0.1 

^•^- 0.1 



Table 2 

Summary of inference results for the strong damping case for the 
three alternative models by employing the TTTB and numerical 
forward solutions. P = 185 s; Td = 200 s; crp = = 30 s. 


TTTB approximation 

Numerical 



M® 

M^- 

M^ M® 

M^ 

M^ 

rAi 

121+2^ 

-20 

120^^® 

^ -20 

121+« I26!|3 

147+39 

-30 

145+27 


5 4+^‘ 
-2.6 

4 0 + 3.4 
- 2.7 

5-3li 8.8)2;5 

c ^+6.8 
-1.8 

10. L« 

HR 

- 0.3 

- 0.2 

0-91:1 0.7^5 

- 0.7 

1.0l;6 


ing cases, Figs.|4^ and|6^, we can see that the impact of the 
strength of the damping on the inversion of this parameter 
is important. Likewise, the probability density functions for 
the transverse inhomogeneity length scale are affected by the 
adopted density model, to a moderate extent under the TTTB 
approximations (see Fig. |6};), but signihcantly for numerical 
forward solutions (see Fig.|6]f). 

Summary data for the inferences performed for the strong 
damping case for the three density models and the analytical 
and numerical forward solutions are presented in Table|2] For 
the TTTB results, we obtain almost the same results when one 
focuses on the median of the probability density function and 
errors given at the 68% credible interval. The situation is dif¬ 
ferent for the case in which numerical forward solutions are 
used. Here, rather distinct results are obtained depending on 
the density model that is employed, although they are masked 
to a great extent when one considers the upper and lower lim¬ 
its of the credible intervals. We can therefore conclude that, 
for observations with strong damping, the inference results 
will signihcantly differ depending on whether TTTB or nu¬ 
merical forward solutions are used. They will also lead to dif¬ 
ferent probability density functions and parameter estimates 
depending on the adopted density model. The differences 
will be more important in the case of the density contrast and 
the transverse inhomogeneity length scale, the two parameters 
that dehne the cross-held density variation. 

Besides the particular properties of the marginal posteri¬ 
ors discussed above, we stress a crucial difference between 
our results and those obtained from th e classic inversion. I n 
the classic inve r sion, a s emph asized bv I Arregui et all ((2003), 
iGoossens et akl (l2008h . and ISoleret alJ (|2()14|) . any point 
along the inversion curve is equally compatible with the ob¬ 
servations. Furthermore, the uncertainty on the measured 
wave properties is not taken into account. This is not true 


under the Bayesian formalism, where the solution is given 
in terms of a degree of plausibility for the different values 
each parameter can take on. Given that observed data values 
can be obtained by different combinations of parameters and, 
considering uncertainty as well, some combinations are more 
probable than others in causing the observed values of period 
and damping time. The calculation related to how many times 
alternative parameter combinations produce a given data real¬ 
ization is contained in the likelihood function. When com¬ 
bined with the prior, the resulting posterior gives the distribu¬ 
tion for the grade of belief on a given parameter value condi¬ 
tional on observed data. 


5. MODEL COMPARISON 

The second level of Bayesian inference is model compari¬ 
son, which enables us to compare the plausibility of the three 
alternative density models to explain observed data in pe¬ 
riod and damping time. Recall that we use the names M', 
with i - S,L,P, corresponding to the sinusoidal, linear, and 
parabolic density prohles used in this study. For each assumed 
prohle, the marginal likelihood. 


p(d\M‘) = 


f 


p{d,e\M‘)de 


f 


p(dlff,M‘)p(fflM')dff, (10) 


is a quantity that provides us with the probability of the ob¬ 
served data d, given that the model M' is true. It tells us how 
well the observed data are predicted by the particular model 
M'. When two competing alternative models, M' and M^, are 
compared one by one to explain a given set of observations, 
the plausibility of one model over the alternative is given by 
their posterior ratio. This can be computed by hrst rewriting 
Bayes’ theorem in terms of the probability of a given model, 
conditional on the observed data as 


p{M\d) oc P{d\M)p(M) (11) 

and applying it to two competing models to compute their 
posterior ratio as 


p{M'\d) ^ p{d\M') p(M') 
p{M^\d) p(d\M^) p(Mj) 

Considering that both models are equally probable a priori, 
p(M‘) - piM-i) - 1 /2, the posterior ratio calculation reduces 
to the evaluation of the ratio of marginal likelihoods for both 
models. The comparison then consists of the computation of 
the so-called Bayes factor. 


BF-J - 

pidiMjy 


(13) 


The magnitude of the Bayes factor is a measure of the rel¬ 
ative plausibility between models. To translate the obtained 
magnitude s into levels of evidence an empirical scale is used 
(see e.g., iJeffrevsl 119611) . We use the empirical table by 
iKass & Raftervl(ll995ll . which assigns evidence for model M' 
against model that is minimal evidence (ME) to values 
of 2 logiBF'y in between 0 and 2, positive evidence (PE) to 
values in between 2 and 6, strong evidence (SE) to values in 
between 6 and 10, and very strong evidence (VSE) to values 
of 2 logiBF'y larger than 10. 

We have hrst computed the marginal likelihood for the three 
alternative models under consideration, by applying Eq. (doll 
to each of them. As an example. Figure |7] shows the result 
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Figure 7. Two-dimensional surface plot representing the marginal likelihood 
for the sinusoidal density model, p{d\M ^), in the observable parameter space 
(P, Td). The magnitude at each point is a measure of the degree of plausibility 
of the model for the observed data. The uncertainty on period and damping 
time is crp = = 30 s. 

for the sinusoidal model, . A given pair of observed pe¬ 
riod and damping time would correspond to a point over the 
displayed two-dimensional surface. The magnitude of the 
marginal likelihood at such a point is a measure of the good¬ 
ness of this particular model in explaining those observed pe¬ 
riod and damping time values. The calculation already takes 
into account the uncertainty on the measured wave proper¬ 
ties. We see that certain combinations of period and damp¬ 
ing time have a larger marginal likelihood than others. This 
means that when performing the integral in Eq. (fTOl) . summing 
up all the possible parameter combinations that could produce 
the observed data, those combinations have a larger plausibil¬ 
ity. The calculation requires the use of that particular model, 
assumed to be true, the selection of particular values for the 
uncertainty on the data, and the integral of the product of the 
likelihood function with the prior over the assumed ranges for 
all parameters; hence, all the available information is used in 
a consistent way. A similar calculation was carried out for the 
remaining linear and parabolic density models, and . 
Once this is done, Bayes factors can be computed considering 
ratios as given by Eq. (fT3l) and the relative plausibility of one 
model against another can be assessed. 

Eigure [8] shows the obtained results. Each panel displays 
contours for the regions in period and damping time space 
where the evidence for one particular model against another 
alternative reaches a given level of evidence. These levels are 
based on the corresponding values of 2 logl BF'Q and the as¬ 
sociate d amount of evidence according to the lKass & Rafter^ 
(119951) table. Regions in white indicate period and damping 
time combinations for which minimal evidence is obtained 
and hence the result of the model comparison is inconclu¬ 
sive. Then, different levels of gray shading point to regions in 
which positive, strong, and very strong evidence is obtained, 
with the level of evidence increasing with the darkening. The 
top row panels correspond to results in which the TTTB ap¬ 
proximations to the forward solutions have been used. The 
bottom row panels correspond to results in which numerical 
forward solutions have been employed. 

Considering first model comparison results under the TTTB 
approximations, Eig. [8^ shows the regions with different lev¬ 
els of evidence for the comparison between the sinusoidal 
and the linear profiles. Note that, by construction, regions 


in the observed parameter space where BF‘^ and BF^‘ reach 
the different levels of evidence are mutually exclusive, since 
\Qg(BF'^) - — \Qg{BF^‘), and cannot overlap. In the com¬ 
parison between the sinusoidal and linear profiles only BF^’’ 
reaches significantly positive values, with the evidence being 
positive, strong, and very strong. This means that the evi¬ 
dence supports the linear density model instead of the sinu¬ 
soidal density model. However, the regions where this occurs 
correspond to combinations with short damping timescale, in 
comparison to the period (see Eig. [8^). A similar conclusion 
can be reached in the comparison between the sinusoidal and 
parabolic density models (see Eig. [Sti) and the comparison 
between the linear and parabolic density models (see Eig.[8j;). 
In the first case, we obtain positive and strong evidence for 
the parabolic profile. In the second case, the evidence sup¬ 
ports the linear density model instead of the parabolic profile, 
with positive, strong, and very strong evidence. In both cases, 
this happens again for combinations of period and damping 
time indicative of very strong damping. In summary, model 
comparison under the TTTB approximations enables us to 
draw conclusions about different levels of evidence among the 
considered models, but this evidence is only appreciable for 
strongly damped oscillations. If that were the case, the linear 
and parabolic profiles are preferred in front of the sinusoidal 
profile. Among them, the plausibility of the linear density 
model is larger than that corresponding to the parabolic den¬ 
sity model. 

Turning now to the case in which numerical solutions for 
the forward model are used, outside the TTTB approxima¬ 
tions, the results are strikingly different (see Eigs. [8}l-f). In 
appearance, Eig. [8^ and d, corresponding to the compari¬ 
son between the sinusoidal and linear density models, look 
rather similar. The difference is that, in the second case, re¬ 
gions with significantly positive values of \og{BF^^) are be¬ 
ing plotted, i.e., those combinations of observed period and 
damping time for which the evidence supports the sinusoidal 
model. The model comparison result using the numerical so¬ 
lutions is therefore contrary to the one obtained under the 
TTTB approximations, where the linear density model was 
preferred. The situation is more involved when we compare 
the sinusoidal and parabolic density models. In this case, see 
Fig. there are regions in observed parameter space for 
which evidence supportive of both models is found. Those 
regions where the evidence supports the sinusoidal density 
model are contoured with solid lines, with the evidence being 
positive, strong, and very strong. The region contoured with 
dashed lines corresponds to period and damping time values 
for which positive evidence for the parabolic density model is 
obtained. Einally, in Eig. [8f, the result from the comparison 
between the linear and parabolic density models is shown. In 
this case, the result is similar to the one obtained under the 
TTTB approximations (compare Eigs. [8}; and f) and the evi¬ 
dence supports the linear profile instead of the parabolic den¬ 
sity model. As before, model comparison using the numerical 
forward solutions enables us to draw conclusions about dif¬ 
ferent levels of evidence among the considered models, but 
this evidence is only appreciable for strongly damped oscilla¬ 
tions. The sinusoidal and linear density models are preferred 
in most of the regions in observed parameter space, with a 
little region where positive evidence for the parabolic density 
model instead of the sinusoidal profile exists. 

6. MODEL AVERAGING 








Arregui, Soler, & Asensio Ramos 


10 








Figure 8. Contour plots of the regions with different levels of evidence according to Bayes factor values from the comparison between alternative density models 
in the two-dimensional observational parameter space {P, rj). ME: minimal evidence; PE: positive evidence; SE: strong evidence; VSE: very strong evidence. 
Top row panels correspond to the analysis using TTTB forward solutions. Bottom row panels correspond to the analysis using numerical forward solutions. The 
density models being compared are: (a) and (d): sinusoidal vs. linear; (b) and (e): sinusoidal vs. parabolic; and (c) and (f): linear vs. parabolic. 


Let us summarize the results gathered so far. In Section |4l 
we have shown that the Bayesian solution to the inference 
problem enables us to obtain well-constrained marginal pos¬ 
teriors for the three parameters of interest. However, the 
obtained results depend on the particularly adopted model 
for the density variation at the nonuniform transitional layer. 
These differences are small and mainly focused on the poste¬ 
riors for the transverse inhomogeneity length scale in the case 
with moderate damping, but can be significant and even affect 
the inference of the Alfven travel time for cases with strong 
damping. The use of numerical solutions for the forward 
problem instead of the analytical estimates under the TTTB 
approximation leads to more marked differences in the cor¬ 
responding inversions. In Section |5] we have shown that the 
application of model comparison techniques leads to different 
levels of evidence for one model to be preferred to another, but 
the regions over observable parameter space where this hap¬ 
pens are located at combinations of period and damping time 
corresponding to very strong damping regimes. This means 
that it might be difficult to perform such assessment with well- 
differentiated degrees of plausibility among the competing 
models, because the evidence for one model to be preferred 
over another is not strong enough for many combinations of 
observed period and damping times. The question then re¬ 
mains as to which one of the inferences is to be preferred. 

The period and damping time values for the moderate and 
strong damping example cases fall in regions where the evi¬ 
dence is minimal in the comparisons shown in Fig. [8] How¬ 
ever, the evidence for each model given a particular pair of 
values for period and damping time is different. For exam¬ 
ple, a calculation of the Bayes factors for the moderate damp¬ 
ing case, using the numerical forward solution and taking the 
sinusoidal model as the reference model, gives - 1; 

BF^ - 0.67; and BF^^ - 1.08. This may lead us to chose 
the parabolic density model as the most plausible one in view 
of data and to therefore adopt the corresponding inference as 
the most reasonable one. However, the differences in Bayes 
factors are not large enough to support this conclusion. 

To solve this problem, a further step can be pursued by con¬ 


sidering a third level of Bayesian inference, namely, model 
averaging. Bayesian model averaging is a procedure to obtain 
parameter constraints that account for the uncertainty about 
the models under consideration. Applied to our particular 
problem, it consists of combining the posteriors obtained for 
each particular density models to obtain a model-averaged 
posterior, weighted with the evidence for each model. 

The model-averaged posterior distribution for a given pa¬ 
rameter 6i, conditional on observed data d and weighted with 
the probability of our set of density models is given by 


p(0i\d)^ ^ p(0i\d,M'‘)p{M'‘\d) 

k = S,LP 

^p(M^\d) Yj BF'‘^p(0i\d,M'‘), (14) 

k = S , L,P 

where in the second equality we have adopted the sinusoidal 
profile, , as the reference model and replaced the mod¬ 
els’ posterior probabilities by the Bayes factors, BF'^^ , with 
respect to the reference model. Obviously, BF^^ = 1. The 
posterior for the reference model, , can be calculated by 
considering that the sum of the probabilities for all three mod¬ 
els must be unity, thus 


p(M^\d) = 


1 

1 + Zk = L.p ' 


(15) 


The obtained result is independent of the model chosen to be 
the reference. 

Figure |9] shows the inference results for the three param¬ 
eters of interest and the two considered cases, for moderate 
(Fig. |9^-c) and strong (Fig. |9}l-f) damping. In all panels the 
model-averaged posteriors using the TTTB and the numeri¬ 
cal forward solutions are compared. In the moderate damping 
case, very similar model-averaged posteriors are obtained for 
the three parameters of interest. Differences are more marked 
in the strong damping case. 

Table [3] displays the summary statistics, with the median 
and errors at 68% credible regions, for each model-averaged 
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Figure 9. Model-averaged marginal posteriors for the three parameters of interest computed using expression <14> . Top and bottom panels correspond to the 
moderate {P = 272 s; rj = 849 s) and strong (P = 185 s; = 200 s) damping cases, respectively, in all cases, crp = = 30 s. In each panel, solid lines 

represent the inference performed using the TTTB forward solutions, while dashed lines are for the inference performed using the numerical forward solutions. 
Summary data for the posteriors with the median and errors at the 68% credible region are shown in Table[3] 


Table 3 

Summary of model averaged inference results for the moderate 
(P = 212 s; Td = 849 s) and strong (P = 185 s; rj = 200 s) 
damping cases. In both cases, crp = = 30 s. 


e, & 

Moderate damping 

Strong damping 


TTTB 

Numerical 

TTTB 

Numerical 

TAils] 

170+24 
i /^_24 

170+11 

^ -19 

141+30 

1 ^ 1-27 



-4.9 

-2.6 


l/R 



0 9+0-5 



1 

1 

1 

1 

BF^ 

0.62 

0.67 

0.71 

1.12 

BF'’^ 

0.89 

1.08 

0.93 

2.42 


inference. Also, the corresponding Bayes factors are shown. 
For the moderate damping case and the TTTB results, the si¬ 
nusoidal model is the one with the largest plausibility, with a 
small difference with respect to the parabolic density model. 
However, when numerical results are employed, the parabolic 
density model has the largest plausibility. Something similar 
occurs for the strong damping case in which the Bayes factor 
for the TTTB result is larger for the sinusoidal model, but the 
parabolic profile turns out to be the most likely when numeri¬ 
cal results are employed. 

The model-averaged posteriors and the summary values of¬ 
fer the most general inference result that can be obtained on 
the unknown parameter values. They take into account all the 
available information, i.e., the prior information, the observed 
data with their uncertainty, the modeling constraints, and the 
evidence for each model in view of data. 

7. SUMMARY AND CONCLUSIONS 

The characteristic timescales and spatial scales for damp¬ 
ing of MHD kink waves by resonant absorption and subse¬ 


quent energy dissipation in coronal waveguides depend on 
the cross-field variation of the mass density. Observed wave 
damping properties can be used to perform the inference of 
the density profile at the nonuniform layer, upon the assump¬ 
tion of particular forms for the density variation in theoret¬ 
ically modeled density tubes. In a recent study. I^ler et alJ 
(|20T3) have pointed out that the particular model adopted for 
the density profile at the nonuniform transitional layer may 
have strong implications for coronal seismology estimates for 
the unknown parameters. The reason is that solving the in¬ 
version problem in the three-dimensional space of unknown 
parameters leads to solution curves that are significantly dif¬ 
ferent. 

In our study, we applied three levels of Bayesian inference 
to the problem of obtaining information on the cross-field 
density structuring in coronal waveguides from observed peri¬ 
ods and damping times of transverse kink oscillations. Three 
alternative density models were considered as a representation 
of the density variation at the nonuniform transitional layer at 
the tube boundary, namely, sinusoidal, linear and parabolic 
profiles. 

The application of parameter inference enables us to deter¬ 
mine the parameters that define the density variation under 
the assumption of a particular model from observed oscilla¬ 
tion properties. The inference result depends on the assumed 
density model. The differences are small when the damping 
of transverse kink oscillations is moderate and when TTTB 
forward solutions are employed. However, significant differ¬ 
ences in the obtained posterior distributions were found when 
the damping is strong and the full numerical forward solutions 
are used. When the posteriors are summarized by consider¬ 
ing, e.g., the median and errors at the 68% credible region, 
the differences are concealed to a great extent. 

By applying model comparison techniques, the marginal 
likelihood for each model and the Bayes factors were com¬ 
puted. They inform us on the likelihood of obtaining given 
observed values for period and damping time and on the rel¬ 
ative plausibility between the considered alternative density 
models in view of data. Our results indicate that although 
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the plausibility for each density model is different and the re¬ 
gions in observable parameter space where the evidence for 
one model against an alternative is stronger can be well differ¬ 
entiated, those regions correspond to period and damping time 
values that imply very fast damping. Therefore, a categorical 
assessment among models, based on positive, strong, or very 
strong evidence, cannot be made for typically observed damp¬ 
ing regimes. Furthermore, as analytical TTTB and numerical 
forward solutions predict different period and damping times, 
the regions in observable parameter space for which the evi¬ 
dence supports a given model are different and lead in general 
to dissimilar conclusions. 

Even if the evidence analysis does not permit to choose a 
particular density model, for any given observation the Bayes 
factors differ. This information on the distinct plausibility of 
each model, in view of data, can be used to perform model 
averaging. The procedure combines the alternative posteri¬ 
ors by weighting them with the individual model evidence. 
Such a calculation makes use of all the available information 
in data and modeling in a fully consistent manner. The result¬ 
ing marginal posteriors are the best inference one can obtain 
with the information at hand. 
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